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1. Introduction. 

Given the solution of a spectral problem, for some differential operator say, the 
calculation of the corresponding functional determinant could be regarded as just a 
computational challenge but there are, of course, uses for such objects. In physics 
they determine the one-loop effective action. In mathematics, for the de Rham 
complex, they occur in the analytic torsion, and elsewhere. 

It is not necessary to have the spectrum explicitly available in order to calculate 
the determinant, but it helps. For this reason many discussions revolve around 
exactly solvable cases and prominent amongst these are the spheres. Some relevant 
brief history was attempted in [1] so nothing more will be said on this, just now. 

The intention of this paper is to present a small contribution to the general 
store of knowledge about spherical determinants, in particular on lens spaces. These 
have played an important part in discussions of analytic torsion, [2]. 

In an earlier work, [3], amongst other things, I calculated the determinants on 
the even homogeneous lens spaces, S 3 /Z g , (see also [4]). In the present work I turn 
to the inhomogeneous case. Curiously the results turn out to look quite different 
when q is an odd prime. A variety of approaches is offered and I also make some 
further technical analysis of the homogeneous case. 

2. Lens space spectrum. 

Thinking of the lens space as the spatial section of an Einstein Universe, I 
consider only the conformal scalar and divergenceless Maxwell vector (coexact 1- 
form) eigenproblems in computational detail. There is a lot of prior mathematical 
work on the spectral problem, as it is relevant for the analytic torsion (the minimal 
scalar is needed here) and the r) invariant, but I shall take the spectral details as 
given in the earlier works, [5-7], since they are in the form which I wish to use. 

The differential operators under consideration are the scalar Laplacian, with 
an addition to make it conformally invariant (in four dimensions), the de Rham 
Laplacian, and, for spin-half, the square of the Dirac operator. 

The fact that the spectrum is composed of squares of integers (up to a scaling) 
means that number theory is almost bound to appear somewhere in the story, and 
this will happen. 

The eigenfunctions are labelled by two angular momenta, L and J where 
J = L+j and j is the spin of the field (j = 0,1/2,1). The massless polariza- 
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tion conditions produce the restriction J = L±j, making L ultimately a sufficient 
label. 

The spectrum on the full sphere of radius a is, 

A+ = 4(.? + n ) 2 ' n=l,2,..., (J = L + j) 

\ (1) 
K = — U - nf, n = 2j + 1, 2j + 2, . . . , (J = L - j) . 

(X 

where n = 2L + 1 e Z and the spectrum has been split into the parts that arise 
from the positive and the negative spectra of the corresponding first order (pseudo)- 
operators. If these parts are to be united, it is necessary to distinguish j = from 
j > 0. The positive and negative eigenfunctions are related by the interchange, 
L <-> J, a parity transformation. The point is that, if j = 0, these modes are 
degenerate and must not be counted twice. Hence it is necessary to introduce a 
degeneracy factor, d(j) = 1/2 ,j = while d(j) = 1 ,j > 0. 

The factoring, S 3 /T, does not alter the positive-negative split and the eigen- 
values are still as in (1) except that the range of n is modified. The degeneracy will 
contain this information. 

The deck group, T, has left and right actions with typical elements, 7 L and 7^. 
The degeneracy takes the SO(4) character form, [7], 



d(L,J) = ± X {L \il)x {J) (1r), (2) 
in terms of the SU(2) characters, 

X (L) (7) = 



sin^L + l^ 



sin 9 1 

where 6* 7 is the 'radial' angular coordinate labelling the group element, 7. 

The freedom allowed by the non-trivial fundamental group to twist the field 
by a representation, Horn (T, U (1)), has not been incorporated here. This would be 
necessary if the torsion were under consideration. 

The spectral data is now combined into a (function on S 3 /r, cf [5], by adding 
the positive and negative parts. (Subtraction would give the rj invariant). 

Trigonometry gives, 

. , s a zs d(j) v-^ 2 v^' 1 , , M , n , 

Cs(s) = —— }^ cos g_ cosa 2^ ^{cosn/3 cos ja - cosnacosj/3) , (3) 

a,P n=j 
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where 

a = e R + e L , p = e R -e L . 

From this point on, the expressions do not apply to spin-half because the 
summation variable has been shifted by j to reach (3). 

Initially, I will proceed without placing the analysis in a wider context, i.e. just 
from a pedestrian calculational viewpoint. 

As in [5] , the (-function is written in terms of the simplest Epstein (-function, 
defined by 

oo 

(*)= J2 \m + g\- s e 2M , (4) 

m= — oo 

with the understanding that, if g — 0, the m = term is omitted. Then, 




a 2s d(j) 



E 



1 



cos P — cos a 



cos ja Z 





P/2n 



(2s) - cos jf3 Z 





a/2n 



(2s) . 



(5) 

A remark of possible interest is that, in our previous work, [3], the ^-function 
on one-sided (homogeneous) lens spaces was reduced to a (^function on a factored 
two-sphere. In the case under study here, the Epstein function, (4), is a twisted 
^-function on the one-sphere, if g = 0. 

I also note that Cs(ty = f° r spin-0 but (3(0) = 1 for spin-1. This will come 
up later. 

The advantage of the Epstein expression is the existence of a functional relation 
that allows (5) to be replaced by the image form, 



_ a 2s d(j) ^ s _ l/2 T(l/2-s) 
^3\ b ) — FFi n 



x 



E 

a,(3 



1 



cos (3 — cos a 



T(s) 
cos jet Z 



/3/27T 





(1 - 2s) - cos jpZ 



a/2n 




(1-2*) . 



(6) 



So far everything has been for a general V. The case of lens spaces, L(q; A 1; A 2 ), 
is covered by the choice of angles 



2tv q ' 2tv q 



(7) 



where p, = 0, . . . ,q — 1 , labels 7. v x and v 2 are integers coprime to q, with A]^ and 
A 2 their mod g inverses. 
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By an appropriate selection of a q-th root of unity, it is possible to set v x = 1, 
i.e. X 1 = 1, without loss of generality. Any pair, (z/ l5 v 2 ), can be reduced to (1, v) 
by multiplying through by the mod q inverse of v x . 

The simple, one-sided lens space, L(q;l, 1), corresponds to setting v = 1 so 
that # L = 0, fi = lirp/q and the (^-function becomes a derivative of an Epstein 
function. This is not the case I am interested in just now. Indeed, if the method 
to be presented later is to run smoothly, it is necessary that the denominator, 
cos (3 — cos a, should never be zero. This puts conditions on q and v which can be 
simply, but not uniquely, satisfied by choosing q to be odd prime, when all v from 
2 to q — 2 are covered. I will do this from now on for organisational convenience. 

The p = value corresponds to the identity element of T, and is best separated. 
For the two spins, the identity ^-functions are, 

Q d (s)=a 2s ±-( R (2s-2), j = 

2 (8) 
Cs d (s) = a 2s ^ {( R (2s - 2) - ( R (2s)) , j = l. 

These differ from the full sphere expressions only by the 1/|T| volume factor. 



3. Lens space determinants. 



On differentiating (5) at s = one encounters the right-hand side bracket 
evaluated at a point where the Z : s have a pole. For spin-0 these cancel, but not for 
spin-1 when they will combine with the l/r(s). There is no problem with this, but 
I proceed in an alternative fashion eliminating the pole using a limiting procedure. 

If g 7^ 0, the Z of (4) has no pole at s = 1. So I insert a non-zero g and let 
g — > near the end after the differentiation with respect to s. It is necessary to 
allow for the extra term at m = introduced in this way. Hence, instead of (6), 



a 



2 s 



d(j)^_ 1/2 r(i/2- s ) ^ i 

r| T(s) ^ cos (3 — cos a 



x lim cos ja Z 



(3/2n 
9 



(1 - 2s) - cos j/3 Z 



a/2n 
9 



(l-2s) 



(9) 



d(j) sr-^ cos j/3 — cosja 



lim - 

, ; „ cos (3 — cos a V g 

a, (3 



2s 
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Now I can employ a formula of Epstein's, valid for v and q integral and coprime, 



[8], 



vjq 
9 



q-l 



(1) = -2 e- 2 ^ k+9)v / q log sin (n(g + k)/q) 



(10) 



k=0 



to give 



lim Z 

3^0 



vjq 
9 



q-l 



(1) = -2 V e - 2ntkl/ / q log sin (irk/q) - 2 lim log sm(7rg/q) , (11) 

' J q^Q 



fc=l 



which is a real quantity and so the exponential can be replaced by a cosine (set 
k — > q — k). 

I note that Epstein's derivation of (10) assumes only the standard summation, 



1) = -21og(2sin7r/i) , 



and also that it can be used to streamline some of Ray's algebra, [2]. 

Since the bracket on the right-hand side of (6) is now finite at s = 0, in 
evaluating the derivative at 0, (^(0), one needs differentiate only the l/r(s) factor, 



d(j) ^ ^ cos ^™ cos ^EEfc _ C os ^£ cos nk 
CM = -2^1 Y 1 & ^ q — log sin 



q *-ifr! cos ^2- cos ^ q 

p=l k—\ q q 

Mi) pi cos 2™ - cos 2™! 

+2^E — Tfe log v + a" (o) 

<Z cos — - — cos — - q 

p=l q q 



(12) 



The second term on the right is the residual effect of the 'zero mode' for spin-one 
and the last term is the contribution of the identity element from (8) . 

The appearance of the radius, a, in the spin-1 expression reflects the non- 
vanishing of £3(0) and the resulting scaling dependence. 

I write (12) cosmetically as, after a few cancellations, 

9-1 irk 9 

C 3 (0) = -J2 ^(°) lo S sin — + -&("2) , 3 = 
k=i 9 9 

C 3 (0) = - J] ^(!) log sin — + -C R (-2) + 2 log(27ra) - J log 2g , j = l 

k=i Q Q Q 

(13) 
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z' 



(0) = - 7 - log27r - V(^) - V(l - ^) , (15) 



where the (rational) coefficients, A fc (j), are defined by 

p=l q q 

I also remark that the same result, (12), follows without a special limit by 
using the remainder at the pole in the Epstein (-function. This is a dilogarithm, 
ip, function, an equivalent statement being, 


h 

which could be used directly in (5). 

One can then employ Gauss' famous formula for tp(p/q), or better, a formula 
that appears during a proof of this relation, [9] p. 146, which I reproduce, (see also 
[10] p.13), 

</>(-) +V>(1 - -) = -2 7 -21og2 + V cos (^) log 2 sin ^ . (16) 
q v q f-^ q k 

k=l 

One could also take the attitude that using the Epstein equation, (10) provides a 
neat proof of (16), and hence of Gauss' formula. 

The functional determinants are conventionally defined by T> (j) = e~^ 3 ^°\ A 
convenient quantity is the ratio R q , to the q-th root of the full sphere determinant, 

V q (0) [( ^ /2] ( . 2 nk^) 

(17) 



u«)=^r q = n 



sin 2 



(2>i(0)) 

w=7^w=(^) n 



V q (l) ( q \ 2A ' " ^ ' ( . 



na , 

' k=i 



sin 2 



where the symmetry A q _ k (j) = A k (j) has been used to fold the product. 

In any specific case these quantities can be computed. The first non-trivial q 
for which the formula applies is q = 5, (v = 2,3). This is because for q = 3, the 
possible value of v = 2 is such that v = 1 mod q and we are back to the one-sided 
(homogeneous or diagonal) case when a different method is needed. The method 
does not work for q = 6 so the next example is q = 7, {v = 2, 3, 4, 5). 

With this in mind, the particular, //-independent values of the A k , 

A 1 (0) = 1--, A 2 (0) = -- , A 3 (0) = 1-- 

q q q (18) 

A 1 (i) = o, A 2 (i) = 2 -^^. 
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can be used to rewrite (17), 

\(<?-5)/<Z/ \(<?-9)/<? 

4 sec 2 2L 1 fe = 4 



for g > 5 and 

for g > 3. The algebraic prefactor is the zero mode effect. 

The above expressions are useful numerically and I treat them, for now, purely 
in this light, as I do equation (14) for the A k : s. The next section has some further 
analysis of the A k quantities and an alternative route to the equations. 

I give a few low values and plot a graph of W = — log_R 2 9(0) against v in the 
spin-0 case, 

— — J ^ 0.46304135, v = 2, 3 

^ 0.3681520, i/=l,4 

R 7 (0) « 0.3212271, ^ = 2,3,4,5 

« 0.1679911, 1/ = 1,6. 
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30 . 



20 . 

W 

10 . 




, , , , , 

' 10 20 ' 30 

twist 

figl. W = — logi?29 f° r conformal scalars on a lens space of order 29 
for twistings v or, equivalently, A, from 1 to 28, (Az^ = 1 mod 29). 

The v = 1 and v = 28 values have been included for completeness. They were 
calculated from the expressions developed in section 5. 

4. Formal elaboration. 

In tune with my general policy of developing alternative, sometimes equivalent, 
approaches, I return to the starting form of the lens space ^-function, (3) ? which, 
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now taking in the definition (14), reads, 



n=l 



n 



2 s 



(21) 



The second term, the non-identity part, is a Dirichlet series. 
The A n 's, are obviously periodic from (14), 

A n+q (j) = A n {j) , 



(22) 



and it is traditional in such cases to break up the sum according to residue-g classes, 
setting n = Nq + k, where N E Z and < k < q. 
Prom (14), 

4)(0) = 0, A (l) = 2/q-2. (23) 

Since A n = A k by (22), this means that k can be arranged to run from 1 to q — 1. 

The non-identity ^-functions, which are the major technical problem, then 
become, in standard fashion, 



^ 1 



q-l 



= (a/q) 



2 s 



(a/q) 



(a/q) 



2 s 



2 s 



k=l N=0 

q-l 

A (j)C R (2s) +J2Mj)C R (2s,k/q) 
k=i 

(q-l)/2 

A (j)( R (2s)+ A k (j)U R (2s,k/q) + ( R (2s,l-k/q) 

k=i \ 

(«-l)/2 



A (j)( R (2s) + J2 MJ)Z 



k=i 



k/q 




(2s) 



(24) 

which constitutes, perhaps, a neater expression than the ones in sections 3 and 4. 
The reflection symmetry, from (14), 



A q - k (j) =A k (j), 



(25) 



has been used in (24), and earlier. 

From general principles, e.g. [11], there are no poles in the non-identity 
function, for fixed point free actions, V. From the behaviour at s = 1/2 of (24) this 
implies the sum rules, 

q-l 

k=0 
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which, together with (23) and ( R (0,w) = 1/2 — w, yield, correctly, the values, Ca(0), 
of the total ^-function as and 1 for j = and j = 1, respectively. 

The local isometry of S 3 and S 3 /r can be further exploited through the small- 
time expansion of the heat-kernel, the coefficients of which are integrals over local 
geometrical invariants. They are thus related simply by a volume 1/|T| factor and 
are determined by the identity (^-functions, (8). The non-identity contributions 
must vanish. (Actually, in the present case, only the first, volume, term in the 
expansion exists anyway.) The coefficients are proportional to £ 3 (— n), and 
so, from (24), setting £g omd (— n) to zero, there follows the sum rules, 

q-l 

J2Mj)B 2n+1 (k/q)=0. (27) 
k=i 

In fact, this is only a check because the expressions on the left vanish identi- 
cally in view of the reflection symmetry, (25), of the A k and that of the Bernoulli 
polynomials, 

B n (l - x) = (-l) n B n (x) . 

There is, however, some valuable tactical information that can be extracted from 
(27). For example, set n = and make use of (26). Then the sum rules reduce to 
the moments, 

q-l q-i ^ 

5>A fc (0)=0, ^2kA k (l) = -(q-l). (28) 

k=0 fc=0 

The fact that the coefficient of l/(n/a) 2s in the (^-function is a degeneracy, 
and hence integral, gives some general information about the A n 's. For spin-0, the 
identity part of the ^-function contributes a (non-periodic) factor of n 2 /q to the 
(total) degeneracy, hence, 

2 

A n (0) = D n — ^— , (D n G Z), or qA n (0) + n 2 = 0, mod q . (29) 

The values (18) fit this pattern. The degeneracy, D n , depends on q and the 
twisting, z/, generally. As a function of n, the image part, A n , oscillates about zero 
and for large eigenvalues D n is dominated by the n 2 term, which is the full sphere 
value divided by q, a volume factor 2 . 

2 The n 2 /q is the Weyl term. Using the sum rule, a sort of discrete eigenvalue counting function 

is N M = Y^n=i D ^ where M = m (l ~ !)> g ivin g N M ~^ M 3 /3q as M -> 00. M 2 /a 2 is the 
eigenvalue, A, and we see that for large M, N M — > X 3 / 2 \A4\/6ti 2 , the Weyl asymptotic law. 
Such a crude argument works only for this leading behaviour 
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As a further check, one can also evaluate the Casimir energy on the lens space 
Einstein Universe, TxS 3 /r, this way. The values agree, as they must, with our ear- 
lier ones, [5,7]. Non essential technical comments about this are therefore relegated 
to Appendix B. 

My main objective in this paper is the derivative, £3(0), and thence the deter- 
minant. This is easy to find by the approach of this section using Lerch's formula, 

C R (0,w) = \og(T(w)/V^), (30) 

as in Ray's derivation of the torsion, [2], or in the evaluation of the class number 
for quadratic forms of positive discriminant. 

The answers are those given earlier, e.g. in (17). In showing this, one needs 
again the sums, (26). 

The main point I wish to make now is that, if q is odd prime, the determinant 
formulae in (17) can be written in terms of units, e fc , of the g-th cyclotomic number 
field defined by 

/ (l-Cj-)(l-O V /2 sin(rt/ g ) 
"~ U -<.)(! -ffV "-»(-/«)' k - 2 - [ "- 11/2 ' (31) 

where ( q = ex.p(2iri/q) is a primitive q-th root. See Hilbert's Zahlbericht, [12], and 
also Borevich and Shafarevich, [13], p. 360. Franz, [14], uses the square of e. 
Using (26) and (18), the results are, 

(9-l)/2 

«.«>)= n "* <0) 

k=2 (32) 

/ . 1 \ 2(1-1/q) («-l)/2 V ' 

V ' k=2 

where qA k (j) G Z. 

The sum rule, (26), plays a simple but vital role in the derivation of these 
formulae. 

Because of the differing powers, A k , these expressions are not invariant under 
all conjugations of the field. 

It is a theorem that the product of units is a unit and also that a rational power 
of a unit is a unit, e.g. [15] §§ 90, 91, 105. Hence, formally, the ratios in (32) are 
cyclotomic units (up to a factor for spin-one). 
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I remark that the corresponding spin-1/2 formula is, 

2q-l 

W2)= II ^ fe(1/2) ' ( 33 ) 

fc=l,3,... 

where e k are 2g-th cyclotomic units, and the A k (l/2) are the (partial) degeneracies 
arising from the non-identity actions. Their form is given in [7] and again we have, 

2q-l 

A(V2) = 0, 

fc=l,3,... 

corresponding to Cs(0) = for massless spin-half. 



5. Return to the homogeneous case. 



In [3], the homogeneous determinants were computed using an expression for 
the lens space ^-function in terms of that on the orbifolded two-sphere. In this 
paper, for variety, I wish to take the spin-0, one-sided Epstein form, [5], 



CM 



a 



2 s 



y— 



d 



2\T\ ^ sin6> 7 86 1 





6J2tt 



(2s) 



(34) 



further to obtain an alternative expression. 

There may be a certain amount of overkill in this, but I have found it helpful 
to have a number of equivalent expressions to hand, if only for numerical peace of 
mind. In view of the many identities and relations for the (^-functions, Bernoulli 
polynomials etc., one should expect several versions of the same thing. 

As before, the image form is, [5] , 



c 3 ( s) = -£>-^ i w 2 - s) 



2in 



V{s) 







sin 9 89^ 

j 7 7 



9J27T 



(1 - 2s) , (35) 



which can be re-expressed in terms of the more familiar Hurwitz ^-function, 



T(s) 



(36) 



This general formula can be used to calculate any quantity required. There 
is, however, a slight awkwardness whenever 9 = 7r, since both numerator and 
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denominator in the summand vanish. This happens in even lens spaces in which 
case this term, as well as the identity one, is treated separately. This can be done 
by extracting the Z 2 C~f unc tion scaled by a volume factor. 
Hence, the formula suitable for even q is, 

( 3 (s)=a 2s - q (l-2 2 - 2s )U2s-2) 

^-3/2 r(3/2 _ 8) i / 

+ — f5T" £ —(^7^) - 2s ' -) - c«(2 - 2 S , - -) 

(37) 

Further manipulations would take us to the expressions in [3,16]. 

Sometimes it is not necessary to make this step. For example, when evaluating 
the Casimir energy this way, by setting s = —1/2 in (36), the Hurwitz functions 
are evaluated at an odd argument and combine to a sum that gives a trigonometric 
result which cancels against the sin 9 on the bottom. The previous answers in terms 
of sums of cosecants, [5,3], can be obtained in this very roundabout, dual fashion. 
I offer some further remarks in Appendix A. 

For completeness, the odd q formula is, 

Cs(*) = yC B (2«-2) 

i j^rgg-jpy i / (2 _ 2s E) (2 _ 2s1 _p_ ) ) 

(38) 

The derivatives at 0, (^(O), follows more or less as before, 
GO) - -^C fi (3) + -L '^—L-^^E) _ Cr(2 , !_£)), (39) 
for odd q and, for even, 

£(0) = 3^C*(3) + ^- £ ■ * , J c R fr-)-C R frl--)), (40) 
7T 5 27rg ^— J sm(2irp/q) \ q q J 

which constitute decent, numerically amenable forms. 
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5. Discussion 



The main results of this short paper are the explicit expressions (17) and (32). 
These yield the determinants on the lens space, L(q; 1, v), with a restriction on the 
coprime q and v which can be met, sufficiently, by restricting q to be odd prime. 

Figure 1 graphically summarises the numerics. It obviously exhibits the lens 
space homeomorphism L(q; 1, A) ~ L(q; 1, A') when A' = ±A (modq) and the home- 
omorphism when A' A = ±1 (mod q) can also be observed, e.g., L(29;l,8) ~ 
L(29;l,ll). These equalities explain the appearance of the graph as a series of 
extrema and therefore account, partly, for the startling similarity with the plot of 
the corresponding Casimir energy in [7]. 

The fact that the image contribution to the determinant is a cyclotomic unit 
is presented simply as an amusing formal identity. 



Appendix A 

I give some further analysis relevant for the homogeneous lens space conformal 
(■-function, another method of obtaining which is to apply Plana summation directly 
to the sum definition, 

°° 77 gjfl 77$ 

C 3 (s, rn) = J^^2 S in# 7 (n 2 + m 2 )* ' 

where, because I can, I have added an extra mass-like term. Basic contour manip- 
ulations, cf [17], lead straightforwardly to the 'real' integral for the non-identity, 
'image', part of the (-function, in the range s < 1, 

Cr -"(, m) = 4a- sin „ £ f , ^ ~ " , < ft, < 2, . 

p Jm \V ' —m z ) s sin v J sinh ny ' 

I now set m = when the integral, [18], [19], 3.524.1, yields (36) which can 
now be extended to all s. Furthermore, [19] 3.524.10-16, enable one to obtain 
trigonometric formulae when s = —1/2, etc. This is yet another longwinded way of 
deriving the cosecant sums in [5,3]. 

Useful general formulae to bear in mind are, 

f~ 2r! inh^ft- 1 ) = _1£ ft 
J yy sinh Try 2 d6 2r 2 K 1 
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JO 



and 

9 „ , i sinhy(# — it) 1 d 2r 9 . 

dyy 2r+1 ^7 ^ = -— ^cosec-. 42 

v yy sinhTry 2 d9 2r 2 v ' 

Dividing by sin#, the right-hand side of (41) is a polynomial in cosec 2 #/2 which 
can be determined by recursion. 3 

Expressed in terms of the series forms, these equations are equivalent to those 
of Eisenstein [21], (see e.g. Hancock [22] p. 32 Ex.5). Defining, 



oo ^ 



(x + n)9 
one has the connections, 



(2,), (^ + l,x) = -l|z 



(2g,x) = Z 
The polynomial referred to above is 

(2g, x) = ^ ^(-l) fc+1 A 23j2fc cosec 2k nx , 
k=i 

with the recurrence relation 

and A 2g 2 g = 1- Eisenstein remarks that the coefficients are related simply to 
Bernoulli numbers. 

Likewise one has the expression 

9 

(2g + l,x) = n 29+1 cos Tix ^(-l) fc+1 A 2ff+li2fc+1 cosec 2fc+1 7nr . 

k=i 

If there are further operations to be performed, such as summing over the an- 
gles, then these polynomials are not necessarily the best way of proceeding. Leaving 
things as series is sometimes more economical, [6,23]. 



3 As a point of historical interest, the quoted integrals in [19] are all taken from the classic 
compilation by Bierens de Haan, [20]. No references are given in this work, the author referring 
for these to vols IV, V and VIII of the Memoires of the Royal Dutch Academy. Only vol IV 
is available to me. The lists in these volumes are, apparently, even more extensive than the 
mammoth [20]! 
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Appendix B 



I elaborate on the computation of the Casimir energy using the (^-function 
form, (21), and make some relational comments on earlier calculations. 
The non-identity contribution is, 



l Cr ud ( _ l/2) 



2a 



(9-l)/2 



4)0X*(-i)+ E ^0-)(c fl (-i^) + Cfl(-i,i-^)) 

k=i \ Q 1 J 



4a 

q_ 

4a 
~2a 



A (J)B 2 + iq j2 MJ) (b 2 & +B 2 (1- h 
k=i v y y 



(9-l)/2 



.fc, 



A (j)S 2 + 2 £ A(J')5 2 (-) 
fc=i 9 

x »» + ( 'g /2 ^0)(^-^ + i) 

fe=i y y 



12 



(43) 

in terms of Bernoulli polynomials. 

For simplicity, I carry only the spin-0 results further by using the vanishing 
moments (28) to give, 



Q k=0 



(44) 



Evaluation of the A k purely numerically from (14) yields, for (43), or (44), 
values in complete agreement with those computed in [7] which were there expressed 
in terms of generalised Dedekind sums. By summing over n first, [5,7], one comes 
quickly to the angle sum form, 



-Q onid (-l/2) 



16qa 



1 9-1 

E cosec 2 np/q cosec 2 Trpv/q , 




(45) 



P =i 



which was determined in [7] by means of a recursion technique as an explicit quartic 
polynomial in q whose coefficients were found as numbers, for any given v. 

Arising from the same ingredients, the angle form (45) must be derivable from 
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the moment form (44). Using the finite Fourier transforms, 



q-l 

^ COS 

fc=0 



2nkpu 



q 



q-l 

0, ^ k cos 

fc=0 



2nkpu 



q 



q 

2 



(46) 



q-l 

y~] k 2 cos 

fc=0 



2nkpv 



q 



- cosec 
2 




q 



on (14), with some trigonometry the desired equivalence can be shown. These 
transforms also confirm the properties of the A^s, (26) and (28). 

At the moment I do not have a means of finding the general form for the partial 
degeneracies, A k , and so the present way of calculating the Casimir energy is no 
better, apart from novelty, than simply performing the sum, (45), numerically for 
given q and v. Furthermore, in the present procedure, q has been assumed prime, 
for convenience. 

For j = 0, the denominator in (14) can be divided into the numerator at the 
expense of introducing a product, 




q-l k-1 



but there seems to be no advantage in this. 
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